A Gaussian Model for the Membrane of Red Blood Cells with Cytoskeletal Defects 
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We study a Gaussian model of the membrane of red blood cells: a "phantom" triangular network 
of springs attached at its vertices to a fluid bilayer with curvature elasticity and tension. We 
calculate its fluctuation spectrum and we discuss the different regimes and non-monotonic features, 
including the precise crossover at the mesh size between the already known limits with two different 
tensions and the renormalisation of the bending rigidity at low wavevectors. We also show that 
the non-diagonal correlations reveal, in "dark field", the cytoskeletal defects. As a first step toward 
a non-invasive defect spectroscopy, the specific case of lacking bonds is studied numerically and 
analytically. 

PACS numbers: 87.16.Gj, 87.16. Dg, 61.72.Hh 

The outer walls of biological cells are fluid bilayer membranes made by the self-assembly of lipid molecules in 
water 0. For about three decades, the physical properties of these two-dimensional fluid surfaces, including elasticity, 
thermal fluctuations, phase behavior, topology, shape polymorphism, etc., have been thoroughly investigated in model 
systems [2-EJIjl In real cells, membranes host a large number of proteins and other inclusions, and they are attached 
to a network of filaments called the cytoskeleton 0. In most cases this network is three-dimensional. In the case of red 
blood cells (RBC's), however, the cytoskeleton is a triangular two-dimensional network of flexible spectrin filaments 
attached to the membrane by proteins located at its vertices 0, 0, • Recent experimental and theoretical studies 
have focused on the shape elasticity H 0, El 0, 113 and fluctuation properties 0, Q, 0, 0, 0] of these composite 
membranes. Up to now, the available theoretical models describing the RBC membrane fluctuations either are mean- 
field or describe the cytoskeleton as a fixed plane producing a checkerboard (square) modulated potential . In 
order to obtain a better description, it is essential to release the assumption of a fixed cytoskeletal plane and to allow 
for coupled fluctuations of the membrane and the cytoskeleton. It is also desirable to take into account the triangular 
nature of the cytoskeleton and to determine the effects of the cytoskeletal defects 0, 0] . 

In this paper, we propose a simple gaussian model describing the composite membrane of RBC's. It allows to 
calculate explicitely the whole fluctuation spectrum, even in the presence of cytoskeletal defects, such as lacking bonds 
and seven-fold or five-fold defects. We model the composite membrane (Fig. ^) as a triangular network of harmonic 
springs attached at its vertices to a fluid membrane described by a curvature Hamiltonian plus surface tension |19| . 
The only coupling between the two subsystems is the constraint that the cytoskeleton vertices lie within the membrane 
surface. As in model polymers [2fij . the gaussian character arises from neglecting the excluded volume between the 




FIG. 1: Sketch of the model red blood cell (RBC) membrane. 
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FIG. 2: Network of N springs (dashed lines) in a M x M unit cell of a cytoskeleton with periodic boundary conditions (dots) 
along the basis vectors ti and ta. Here, M — 3 and TV = 26. The extremities of the nth spring have projections R n and H n 
onto the reference plane. A bond is lacking, in order to illustrate a simple situation with a density of defects. 



membrane and the cytoskeleton (thereby treating the latter as "phantom" |23). Depending on the experimental 
conditions, RBC fluctuations may be of thermal origin, or they may be driven by ATP consumption |22t l23|. Here, 
we consider only the thermal, equilibrium situation. We derive the height fluctuation spectrum (|/i q | 2 ) of our model 
system, and we discuss the crossover between the different regimes and the non-monotonic features arising from the 
membrane-cytoskeleton coupling. Since the cytoskeleton and its defects break the translational symmetry, we study 
the extinction conditions for the non-diagonal spectrum, i.e., (hqhk) for k ^= — q. We demonstrate that the latter 
reveals, in "dark field" , the cytoskeletal defects. Detail studies should therefore allow to determine the nature and 
the density of the defects. Here, we discuss preliminary results concerning the spectral signature of lacking spectrin 
bonds. 

Because of the triangular symmetry of the cytoskeleton, we introduce a planar reference frame with unitary basis 
vectors ti and t 2 making an angle of 27r/3 (Fig. |2J). We define the dual basis (t x ,t 2 ) by t Q • t' 3 = 6@, where 6@ is 
the Kronecker symbol. The metric tensor g a p = t a ■ tp has diagonal (resp. off-diagonal) elements equal to 1 (resp. 
— |) and determinant g — det(g a p) — |. Vectors in the plane will be represented by their coordinates according to 
x = x a t a = x a t a (Einstein summation assumed throughout), where x a = g af3 xp with g a ^ = g~p. We describe the 
membrane by its elevation /i(x) = h(x a ) above a reference plane parallel to (ti, t 2 ) and we assume periodic boundary 
conditions such that h{x 1 ,x 2 ) = h(x x + Li,x 2 ) = h(x 1 ,x 2 + L 2 ), Vx Q . Our model Hamiltonian (per unit cell) is the 
following: 

L l n 

h= (I vS^W 2 [f (d a d«hf + Ud«h){d ot h)\+^^[h{K)-KK)) 2 . (i) 

■'° ^° n=l 

The first term, Tto{h], is the standard membrane Helfrich Hamiltonian |l9j| expressed in our non-orthonormal basis, 
with k the membrane bending rigidity and a an externally imposed surface tension. The second term, V[h], modeling 
the elasticity of the spectrin network, represents the stretching energy associated with the N springs that lie within 
the repeat period. In this expression, and are the projections onto the reference plane of the extremities of 
the nth spring (Fig|2J|. We take for R" and R™ the equilibrium positions that they would occupy in the case the 
membrane were flat, since in the limit of small Gaussian fluctuations the in-plane and the out-of-plane fluctuations 
of the vertices are decoupled. We allow for any kind of defect or heterogeneity, therefore we do not assume that 
R" and R" make a regular triangular network. The constant fi n is related to the rigidity k and free-length I of the 
springs (which we assume to be identical biological elements) and to the geometry in the following way: to second 
order in h, the spring elongation energy ifc[{(R n — R rl ) 2 + [h(R n ) — ^(R,,)] 2 } 1 / 2 — £] 2 gives the in-plane contribution 
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|fc(|R„ - R„| - £) 2 plus the term V[h] of Eq. (JTJ, with 



M „ = fc(l-£|R„-R n |- 1 ). (2) 

Because the network may have been stretched or compressed, either mechanically or osmotically, or because of the 
possible defects, the \i n may have either sign. 

Let us define a Fourier transform by h(x x , x 2 ) — {yfgL\Li)~ x l 2 X] q e lqaX where the wavevectors are quantified 
as q a = n a 27r/L a , where n a G 7L. In order to work with dimensionless quantities, we normalize the energies by 
k^T = the in-plane lengths by £ CT = yj n/o and the membrane height h(x a ) by 1/ ' \/fto . Hence ft, q is normalized 



by */ n/ (3 a 1 and fj, n by <r. Unless otherwise specified, all quantities are now dimensionless. We thus have TLo = 
S q k i^q'-'q k^k where G q .k = 5 q+ k/(q 2 +q 4 ) and q 2 = q a q a , yielding the correlation function of the bare membrane: 

For an infinite membrane, this gives a Bessel function of the second kind: G(x) = — iYo(|x|). Because TL is not 
diagonal in q-space due to the dependence in R n and R„, it cannot be directly inverted to yield the full correlation 
function G(x) = (h(O)h(x)). To proceed, we add an external field, TL — > TL + X) q ^q/-qi an d we calculate (/i q /ik) = 
— <9 2 F/(<9/_ q (9/_k)|/=o, where F = — ln[J T>[h] exp(— TL)] is the free-energy. In order to integrate the membrane 
degrees of freedom, we rewrite exp(— V) with the help of auxialiary fields (one for each spring) as 

Ij{d4n cx P ]T (-^ + i<t> n [HK) - KK)]) ■ (4) 

Sums and products over n (and further on rn) implicitely run from 1 to N. Performing the gaussian integral over h 
yields exp(— F) — JY\ n d(j) n exp(—TL), where 

* = E $- \ E(^-/- q ) G ^ ( 5 -k-/-k) , (5) 

n Z q,k 

in which Sq = (^/gLiL 2 )~^ J2 n i4 , n{e~ lq °' R ™ — e~ l9 °^"). The integral over the 4> n is again Gaussian, yielding F[f] = 

-|Eq,k/-qG ; q> k/-k + ^[/],with 



w] = ^E(E c q n) /- q )^(E 

n,m \ q / \ k 



<£%k , (6) 



where 



C («) = - - - (71 

q (V5ilL 2 )i q 2 + q 4 {,) 

Bnm = — S nm + G(R n — R m ) — G(R n — R m ) — G(R n — R m ) + G(R„ — R m ) . (8) 

It follows that the full correlation function r q k = (hqh^) = G q k + Ar q k is given, with no approximation, by 

A N 

r . - q+ k V r< n > Z?" 1 r (m) (91 

q ' ~~ a 2 4- a 4 c q ™ m c k • W 

M 4 n,m=l 

For weak cytoskeletal strains, i.e., ^„ = O(^) <C 1, we have B nni — f-^n^nm 

+ 0(n 2 ), and thus 

Ar - k " -^L^E^ (q2 + q4)(k2+k4) (10) 

where ^„ jq = e _^q ' R '' , — e - i q R ™ . The fluctuation spectrum can thus be obtained analytically for weak strains even in 
presence of defects, or, for larger strains, numerically by simply inverting a N x N matrix. 

Uniform cytoskeleton. - In the case of a uniform dilation or compression in the absence of defects, the points 
R n and R„ form a regular triangular lattice with fi n = fx, Vn, corresponding to a uniform meshsize £ such that 
N = 3LiL2/£ 2 ■ Thus, in dimensional units and for small values of n/a, Eq. Ijlt)|l gives for k = q the following 
fluctuation spectrum: 
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FIG. 3: Fluctuation spectrum of the model RBC membrane with a uniform triangular cytoskeleton, for a tension corresponding 
to \J~a~fk ~ 4.5 £ _1 and a direction of q corresponding to — 0. The case > (resp. ji < 0) corresponds to a stretched (resp. 
compressed) cytoskeleton. We chose \fi\/(J = 0.5 to better emphasize the features of the curves, although Eq. pip is valid in 
the limit fi <C a. All quantities are in dimensional units. 



with 9 the angle between ti and q. For g<( 1 , we get k^T / (|/i q | 2 ) — a e B q 2 + K cS q 4 , with 

a cS ~ CT + ^=A* = v + ^(l-£/0, (12) 

KcS ~ K ~ 32^^ = K+ rl fce(C ~^ ' (13) 

and for q ^> Eq. (flip yields the bare membrane spectrum ksT / (|/i q | 2 ) — a q 2 + nq i . We thus recover the effective 
tension jump a — a e e discussed in Refs. |l4l Il5j . In addition, we obtain the renormalization of the bending rigidity 
K ff due to the cytoskeleton. Taking k ~ 10 -5 J/m 2 ^2] an d £ ~ ^ — 10~ 7 m 24| . we obtain re e g — /c ~ 10~ 20 J, wich is 
comparable but smaller than the usual membrane bending rigidity, in agreement with the measurements of Ref . |25j | . 

The full fluctuation spectrum for large tensions, i.e., a > k/£ 2 , is shown in Fig. 13 There are three regimes: i) a flat, 
tension-dominated regime for q < revealing cr e ff, ii) a crossover region < q < \J a / 1 n, showing the transition 
from (Tcff to a and thus exhibiting a negative slope for large enough fi > 0, iii) a terminal linear regime dominated 
by the bending rigidity, revealing k. In the case of weak tensions, i.e., a < k/£ 2 (not shown), the transition to the 
terminal regime occurs before q reaches £ , hence the change in tension becomes invisible. Since n c g ~ k, one does 
not see either the slope change occuring around q ~ The curves thus look like the dashed ones in Fig. El except 
that the intersection of the asymptotes occurs in the region q < 

Cytoskeleton with defects. - Let us discuss the general conditions under which the couplings (hqhk) are extinct or 
not. Our Hamiltonian can be expressed as TL = J d 2 x d 2 r h(x) H(x, r) h(x + r) = jcPqcPk /i q -ff q ,k ^k- In the absence 
of a cytoskeleton, the translational invariance of the system yields the well-known non-extinction condition q = — k. 
In the presence of a cytoskeleton, this invariance breaks down. For a perfectly uniform cytoskeleton, however, H (x, r) 
remains invariant for discrete translations of the form x a — > x a +£ and the non-extinction condition becomes q a = —k a 
modulo 2tt /£. This weaker condition is still strong, because 27r/£ is very large (macroscopically) . Now, in the presence 
of random structural defects, the translational symmetry vanishes and all the couplings (/i q /ik) appear |2fi| |. 

There are mainly two types of defects: lacking bonds and dislocations made by pairs of five-fold and seven-fold 
defects 0, . A s yste matic study of the signature of these defects according to their density and type is outside the 
scope of this paper |27j |. Here, we limit ourselves to the simplest situation, i.e., lacking bonds, as in Fig.0 We show 
in Fig.0]the corresponding non-diagonal spectrum (i.e., (/i q /ik) with q ^ — k) for q = <7 "t 1 / 1 1 "t 1 1 1 and k = — gt 2 /||t 2 ||, 
as obtained from Eq. (|10[1 after numerically minimizing the in-plane elastic energy in order to determine the positions 
R" and R" of the springs ends. A satisfying analytical approximation of this spectrum (Fig. |2J can be obtained by 
neglecting the lattice distorsion, i.e., by using the positions R" and R% of the regular lattice while setting one of the 
fi n to zero. In this case, Eq. (fTTTft gives |(/i q /ik)| = fi ksTg~2 L~ 2 \sm(q ■ r) sin(k • r)\(aq 2 + «g 4 )~ 2 , with r = 
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FIG. 4: Non-diagonal fluctuation spectrum (/i q /ik) in the presence of lacking bonds, for q ^ — k but ||q|| = ||k|| = q, as detailed 
in the text. The squares (resp. triangles) correspond to one missing bond per unit cell, as in Fig. [5] but for M = 80 (resp. 
M — 20). The values of q are quantified by the periodic boundary conditions. The filled symbols are numerically calculated 
from Eq. IIQU . while the open symbols correspond to the analytical approximation discussed in the text. 



where u is the unit vector parallel to the missing bond. In the q — > limit, we obtain 

(T q 

which shows that the non-diagonal spectrum is proportional to the defect density p = l/(3Af 2 ). 

In conclusion, in the light of this model experimental measurements of RBC equilibrium fluctuations under strain 
could give new informations concerning the structure and elasticity of the cytoskeleton. Since the defects and the 
plasticity of RBC's are believed to be controled by the ATPase activity [l^, it should also be interesting to investigate 
the (equilibrium) non-diagonal spectrum at various stages of the erythrocyte life or in particular conditions, such as 
capillary transport. Finally, it should be possible to study in a Langevin picture the dynamical properties of this 
model system, thereby determining the consequence of non-equilibrium ATP-driven noise on the fluctuation spectrum. 
This would be extremely interesting, as recent experiments reveal strong non-e quilibrium effects [ill El . 
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